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A  SOLVmC  ALGORITHM  TOR  SYSTEMS  OF 
NONLINEAR  ALGEBRAIC  EQUATIONS 

1.  INTRODUCTION 

The  puriK>se  of  this  report  Is  to  draw  the  attention  of  prospective  users  to 
the  possiblity  of  solving  exactly  systems  of  algebraic  equations  that  occur  in 
various  contexts. 

The  first  system  to  be  considered  consists  cf  one  quadratic  equation  in 

n  unknowns  (x  ,  x  ,  .  .  . ,  x  )  and  n-1  linear  equations  in  the  same  variables: 

X  Z  n 


n 


Though  the  following  results  are  valid  (except  where  noted  otherwise)  when  q, 
a,  and  b  are  complex  numbers,  the  results  are  most  useful  for  the  real  case. 

The  system  of  equations  (1)  can  be  interpreted  in  several  ways.  First, 
one  can  look  upon  it  as  arising  from  the  purely  geom'^tric  problem  of  aeter- 
miniiig  tb»  set  of  points  of  simultaneous  intersection  of  an  n-dimensional  quadric 
surlace  with  n-1  hyperplanes.  It  will  be  seen  later  that,  for  arbitrary  coefficients, 
a  solution  does  not  always  exist.  Moreover,  because  of  the  nonlinearity  of  the 
system,  we  can  expect  the  existence  of  more  than  one  solution.  This  is  clearly 
evident  for  the  case  n  =  2  involving  the  real  points  of  intersection  of  a  conic 
section  (ellipse,  parabola,  etc. )  with  a  single  straight  line. 

Second,  a  system  of  type  (1)  with  three  unknowns  arises  in  a  famous 
problem  of  classical  mechanics  (discussed  by  Euler  and  bearing  his  name) 
concerning  the  motion  of  a  rigid  body  about  a  fixed  point.  Golubev^  presents 
a  solution  for  this  case  by  a  method  attributed  to  Hess^  and  later  modified  by 
Shiff.3 
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Third,  still  aDCther  geometrical  interpretation  of  system  (1)  is  possible. 
Assume  that,  in  an  n-dlmensional  space  endowed  with  an  Inner  product  of  the 
form 


n 


it  is  required  to  find  a  vector  x  =  (x^,  .  .  . ,  x^)  of  given  norm  (re:  length)  b 
and  given  projections  b^^  on  n-1  fixed  vectors  Cjjj  =  (Cj^^,  .  .  . ,  Cmnl.  where 
m  =  1,  .  .  . ,  n-1.  The  resulting  system  of  equations  is  then  precisely  that  in 
(1),  where  the  vectors  ,  a,,,^)  are  related  to  the  vectors  Cm  by 


(3) 


in  which  Q  represents  the  matrix 
product  (2). 


)  appearing  in  the  definition  of  the  inner 


4 

Finally,  system  (1)  might  be  relevant  to  problems  of  nonlinear  programming 
where  the  quadratic  form  will  represent  a  cost,  or  objective,  function  while  the 
linear  equations  will  represent  constraints,  or  side  conditions. 


In  section  11  we  describe  the  method  of  solution  for  system  (1)  and  discuss 
its  possible  solutions  in  detail.  In  section  III  applications  of  this  algorithm  to 
se^'Ciral  simple  examples  are  presented.  In  section  IV  extensions  of  the  method 
to  various  systems  —  some  of  which  consist  of  equations  with  nonlinearities  of 
higher  degree  (possibly  trarscendental)  —  are  presented,  and  the  associated 
problem  of  extremizing  a  given  function  subject  to  linear  and  nonlinear  con¬ 
straints  Is  discussed.  In  the  final  section,  the  solving  algorithm  is  investigated 
from  the  viewpoint  of  implementation,  and  arguments  concerning  uniqueness  of 
solutions  are  presented.  In  addltou,  it  is  shown  that  the  solving  algorithm  is 
relevant  to  the  nonlinear  programming  problem. 


II.  METHOD  OF  SOLUTION 

The  solutions  presented  below  are  obtained  by  a  straightforward  method, 
in  marked  contrast  to  that  described  by  Golubev,  which  seems  needlessly 
conaplicated  by  the  use  of  Grammians  and  partial  diSerentiation. 
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We  staxt  by  writing  the  equations  of  system  (1)  in  more  compact  form: 


xQx  =  b 

QfyyX  ~  j  tn  =  i  j  j  n~i  , 


where  the  tilde  ~  denotes  the  transpose  (hence  x  is  a  column  vector). 

To  go  on,  we  assume  that  the  vectors  a^,  m  -  1,  .  .  .  ,  n-1,  are  linearly 
independent.  The  opposite  case  will  be  discussed  later.  This  linear  independ¬ 
ence  of  the  implies  that  there  exists  at  least  one  minor  determinant  based 
on  n-1  columns  of  the  rectangular  matrix  m=l,  ....  n-1;  k=l,  .  .  . ,  n) 

that  does  not  vanish  identically.  This  being  the  case,  we  may  assume,  without 
loss  of  generality,  that  the  variables  Xj^  have  been  already  lra>eled  in  such  a 
way  that  the  first  n-1  columns  of  the  rectangular  matrix  yield  a  nonzero  mino’  . 
Similarly,  no  generality  is  lost  if  the  matrix  Q  is  assumed  to  be  symmetric, 
since  this  can  always  be  accomplished  by  a  suitable  redefinition  of  the  coeffi¬ 
cients. 


In  view  of  these  assumptions,  let  A  denote  the  (n-1)  x  (n-1)  matrix  con¬ 
sisting  of  the  first  n-1  rows  and  first  (n-1)  columns  of  the  rectaxgular  matrix. 
Then  system  (4)  can  be  put  in  the  form 


X  Q.A  =  b 


(5) 


where  x*  =  (x  ,  . 

Vl.n>' 


X  );  b*  =(b  ,  .  .  .,  b  ):aQda*=(a 
n-1  1  n-1  1,  n 


2,  n’ 


By  ’.'irtue  of  c  ir  assumption  of  linear  independence,  the  matrix  A  is 
invertible,  so  that  it  is  possible  to  solve  for  x*  from  the  linear  equations  in 
(5): 


By  traiupof  itlon,  we  have  also 


A*  =  ( tr . 

If  we  partition  the  quadratic  form  to  conform  to  the  definition  of  x*.  we 
can  write  for  the  first  equation  ii  (5): 


X'^CTP  -r  =L.m 

’  'nn 

where  Q*  =  (qj.;  1  =  1 . n-1;  j  =  1,  .  .  . ,  n-1)  and  q*  =  (q^ ,  q2jj,  .  .  . , 

‘In-W-  Subati  tution  of  expressions  (6)  and  (7)  for  x*  and  x*  Into  (8)  yieldn 


a^)Xn 

(9) 


In  both  (8)  and  (9),  the  assumed  symmetry  of  Q  has  been  used.  Moreover, 

since  Q*  is  also  symmetric,  (9)  can  be  rewritten  as  the  following  quadratic 

in  X  : 
n 

Q>C  -i-  ^  =  O  .  (10) 

where 

o<  =  f  <ii) 

fhn  ^ 
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r  =  i-M”  ^ . 

(•3) 

The  two  possible  aolutioos  of  (10)  are  then 

(14( 

From  (14),  we  are  led  to  the  followi'ig  two  [jossible  cases: 

1.  aii  0 

At  most,  two  distinct  solutions  for  Xn  arc  'xj.ssible.  1  hts  implies  the  ex¬ 
istence  of,  at  most,  two  distinct  solutions  x*for  the  original  problem  (as 
dictated  by  (6)  or  (7)).  Exactly  two  distinct  solutions  will  be  obtained  if  and  only 
if  -  ay/  0.  In  the  opposite  case  of  P”  -oV  o,  there  is  onl>  one  solution. 

2.  a  ^  0 

From  (10),  we  obtain 


-f-  ^  —  O  (15 

This  implies  that  a  unique  solution  will  be  obtained  in  this  case  if  and  only  if 

fi/  0. 


The  opposite  case,  p  0,  leads  to  two  possibilities: 

a.  For  >  =  0,  it  is  clear  that  any  x^j  will  satisfy  the  original 

equation  (10).  This,  in  turn,  implies  that  there  are  an  infinite 
number  of  solutions  depending  on  one  parameter  x^  (or  that 
there  is  a  «  -^—manifold  of  solutions).  The  geometric  inter¬ 
pretation  of  this  result  is  that  the  quadric  surface  and  the 
hypierplanes  intersect  in  a  single  "straight”  line. 
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b.  For  y  0,  (15)  will  be  Incoosistent.  Hence,  In  this  subcase  no 
solution  exists  for  x^,  and  nautatts  mutandis,  no  solution  is 
possible  for  the  original  problem.  This  result  raises  no  ques¬ 
tion  if  the  rectangula.'  (n-l>xn  matrix  of  the  a-coe£ficients  has 
only  one  nonvanishing  <n-l)x(n-l)  -  minor.  On  the  other  hand, 
if  there  are  other  nonzero  minors,  one  can  ask  whether  all  of 
them  will  yield  the  same  negative  result.  Though  it  is  difficult 
to  prove  in  general  that  this  is  indeed  the  case,  we  expect,  on 
intuitive  grounds,  that  the  various  nonzero  minors  will  produce 
identical  results. 

From  the  discussion  presented  above,  it  is  clear  that  the  existence  of  a 
real  solution  (when  all  the  entities  <2ntering  the  original  system  are  real)  re¬ 
quires  satisfaction  of  the  condition 

^  O  .  (16) 

In  the  opposite  case,  a  real  solution  does  not  exist. 


IIL  .APPLICATIONS 

In  this  section,  the  algorithm  presented  above  is  applied  to  several  sim¬ 
ple  examples  to  obtain  some  well  known  results. 

A  vector  x  =  (x^,  ....  x^)  is  to  be  found  whose  norm  (length)  is  I  .  In 
addition,  x  must  be  orthogonal  to  n-i  fixed  real  vectors  aQj,  m  =  1,  .  .  . ,  n-L. 
These  constraints  lmpl>  chat,  in  the  terminology  of  the  previous  sections, 

js  ^  J  tn  im  i  .  ^  tt —  1 

D  —  ^  /  *-'nry  ^ 

For  the  norm,  we  assume  the  Euclidean  length 

n 

(18) 

This  implies  that  Q  and  Q*  are  square  unit  matrices  of  n  and  n-l  dimensions, 
respectively,  while  ~  ^  q*  =  0, 
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Substitution  of  these  values  in  (11),  (12),  and  1 13)  yields 

GC  = 

(19) 

P  = 

o 

(20) 

t  = 

- 

(21) 

The  form  of  a  reveals  that  it  cannot  vanish  since,  to  vanish,  the  quadratic 
form  In  a*  of  (19)  must  be  equal  to  -1,  which  is  impossible  because  the  deter¬ 
minant  of  this  quadratic  is  seen  by  inspection  to  be  positive.  This  means  that 
the  quadratic  form  has  a  single  minimum  at  a^  =  (0,  .  .  . ,  0)  and  is  positive  at 
all  other  points.  By  assumption,  a*  ^  0,  the  form  is  positive  definite,  and 
a  >0. 


Use  of  these  resudts  gives  the  solutions 


The  existence  of  these  two  distinct  real  solutions,  as  pointed  out  above,  is  a 
basic  feature  of  the  solution.  The  following  special  cases  of  these  results  illus  ¬ 
trate  this  multiplicity  of  solutions. 

First,  we  specialize  to  the  2-dimen8ional  space,  the  fixed  vector  being 
a  =  (Ui,  U2).  Then,  assuming  Uj^  /  0,  we  have 

j;  /»-'=  -t/u, 

oC  —  i  ^  . 

Then 

=  ±  -i 


(23) 


(24) 
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The  pictorial  representation  of  this  result  is  shown  in  figure  1.  The  solu¬ 
tion  vectors  are  of  length  I  and  are  perpendicular  to  the  fixed  vector  a. 


Figure  1 


Next,  we  specialize  to  the  3-dimensional  case  where  there  are  two  fixed 
vectors  a  =  (u  ,  u  ,  u  )  and  a  =(v  ,  v  ,  v  ),  Assuming  that 


(25) 


we  get 


U,  ti 


v;  V, 


(26) 


and 


+  (t-u,  -  u,  v-Xl^  ■ 


(27) 
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The  two  solutioiis  are  then 

Xt  =  -t 

i -  IT  /ST  . 

Let  choose  t  =A^/o.  The  solutions  of  (28)  are  now 

x  +  =  ^ j  tX '•7-w, Cu.iT^-u^vj) 


(29) 


Wt  .ecogn  Zb,  in  the  expression  for  the  vector  x^,  the  usual  definition  of 

the  cross  product  of  t\  n  vectors  a  ,  a  given  by  the  symbolic  formula 

X  £t 


3,  X  = 


M, 


(30) 


where  i,  J,  k  are  the  unit  vectors  in  a  Cartesian  system  of  coordinates. 

The  geometric  representation  of  this  case  is  shown  in  figure  2.  The  two 

Solutions  are  perpendicular  to  the  plane  defined  by  the  two  vectors  a  and  a  . 

X  M 

It  stK>uld  be  noted  that  the  choice  of  x^  as  the  cross  product  of  Sj^  and  sl^ 
correspor  ds  to  a  right-handed  system  of  coordinates,  while  x~  corresponds  to 
a  left-handed  system. 


9 


1 


TR  5353 


Figure  2 


We  can  apply  the  previous  formalism  to  the  problem  of  finding  the  extrema 
(if  existent)  of  the  quadratic  foimi  xQx  subject  to  the  constraints 


a^x  =  h. 


m  =  - 


V 


n-1 


(31) 


We  substitute  the  results  of  section  II  (with  the  same  msumptions)  given  in  (6) 
and  (7)  in  the  function 


(32) 


which  then  becomes  a  function  of  the  single  variable  x  ,  given  by 

n 

FU„)=  OiX^+ZPX„-i-r.  (33) 
where  a  and  pare  defined  in  (11)  and  (12),  respectively,  and  Y  is  defined  by 


>r= 


(34) 
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A  discussion  similar  to  that  following  (14)  yields  the  two  results  below: 

1.  0 

The  extremum  ox  Is  at  the  point  given  by  f’(x^)  =  0: 


X„  =  —  l^/o^  -  (35) 

This  solution  will  be  finite  if  which  is  the  case  considered  here.  This  value 

of  Xq  will  yield  a  minimum  if  or  >  0  and  a  maximum  if  a<  0.  Therefore,  the 
unique  extremum  for  flxj,  X2,  .  .  . ,  Xjj)  in  th  ?  present  case  will  occur  at 

x=(' 

where 

(  f^/<x)Q^J/\  ^  .  (37) 


2.  g=  0 

In  this  cafe,,  F(x^)  takes  the  simpler  form 

=  2_^x„  -tr. 

Ji  0  ^y,  then  (38)  is  the  equation  of  a  straight  line,  and  it  is  obvious 
that  no  extremum  exists  in  such  a  case. 

If  P=  0,  the  function  F(x^)  has  the  constant  value  Y  throughout  the  region 
of  space  described  by  the  side  conditons  (31),  and  again  an  extremum  does  not 
exist. 
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IV.  GENERALlZATi  JNS 


In  this  section,  we  present  generalizations  of  the  previous  results  in 
several  directions. 


Linear  Independence  of  Vectors  a^Q 


In  the  first  instance,  we  retain  the  assumption  of  linear  independence  of 

the  vectors  a  introduced  in  section  II. 
m 

1.  We  start  with  an  extension  of  our  results  to  the  following  system; 

-i-  X  Qx  =  L 

In  (39),  p  =  (Pj^,  p^,  .  .  . ,  p^)  whore  again  the  p|^  values  are  complex  numbers. 

Under  the  assumption  of  linear  independence  of  the  vectors,  we  can 
use  (6)  and  (7)  in  the  first  equation  of  system  (39).  The  resulting  equation  is 
a  quadratic  in  x  ; 


-hi'  =  O  ,  (40) 

where  a  ,  p  ,  and  V  are  defined  below  as 

-2  T a*A^ ct A' (41, 

'nn  • 

!r-  Q* A" e+  f>*A-'W  -  h  (43, 


12 


TR  5353 


and 


The  existence  and  uniqueness  of  solutions  to  the  system  (39)  follow  the 
pattein  established  for  system  (1)  and  hence  need  not  be  repeated.  As  a  con¬ 
cluding  note  to  this  extensi^.'x.  w&  remark  that  one  might  be  tempted  to  try  a 
reduction  of  system  (39)  to  system  (I)  by  a  transformation  of  coordinates  that 
will  eliminate  the  linear  term  px.  This  Indeed  can  be  accomplished,  but  omy 
In  the  special  case  where  Q  is  Invertible  (i.  e. ,  |  Q  |  ^  0).  For  this  special 
case,  the  appropriate  transformation  is  the  translation 

X  =  y  -h  ^ 

where 


It  should  be  clear  that  the  approach  of  (45)  and  (46)  is  more  restrictive. 

2.  A  second  generalization  is  concerned  with  the  system 

f(Xt , .  =  b 

Useof  (7)  in  the  first  equation  of  system  (47)  yields  an  equation  in  one  variable 
X  : 

n 

FYXn)  = 

The  solutions  of  system  (47)  will  depend  on 
In  particular,  if  the  function  fiy^,  .  . 

■f'  ^ 


(48) 

'he  nature  of  the  solutions  of  (48). 

. ,  X  )  is  of  the  form 
n 

xQx)^  (49) 
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In  which  c  is  a  constant  scalar  and  p,  Q,  x  have  the  same  significance  as  in 

(39),  then  the  equation  for  x  will  be 

n 

L.  /50) 

where  a,  p  ,  and  "V  can  be  easily  written  in  terms  of  e,  p,  Q,  A,  a*,  and  b*. 
Moreover,  if  g(z)  has  an  inverse,  then  the  soi.utiors  of  (30)  will  coincide  with 
the  solutions  of 

In  the  evaluation  of  g  ^(b),  attention  must  be  paid  to  the  branch  chosen  tor  g 
since  the  function  g~-‘^(z)  can  be  multivalued.  Thus,  we  see  that  it  is  no  more 
difficult  to  solve  a  system  such  as  (49>  than  to  solve  system  (39).  This  result 
can  be  of  Importance  in  problems  for  which  the  function  ((x^,  .  .  . ,  Xj^)  may  be 
approximated  by  a  function  of  type  (49).  A  truncated  Taylor  series  for  f  is  a 
special  case  of  this  form. 

3.  A  third  extension  deals  with  the  system 

^  Xn) 

. . 

X  J  rYt  ~  j  t}—2 


where  P,.  Pg-  and  b  are  complex  numbers  and  the  vectors  a^j  are  again 
assumea  to  be  linearly  independent.  This  assumption  implies  that  there  exists 
at  least  one  nonzero  (n-2)x(n-2)  minor  of  the  a-coefn dents.  Again  assuming  a 
labeling  such  that  the  first  n-2  columns  yield  this  minor,  we  can  write 

K*  =  > 
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(54) 

and  A  Is  the  matrix  of  the  nonzero  (n-2)  x  <u-2)  minor. 

Substitution  of  (53)  into  the  first  two  equations  of  system  (52)  results  in 

the  two  simultaneous  equations  for  x  x  : 

n-1’  n 

ir(x„  ,x„)  =(’ 

If  =jft  ■  (55) 

The  points  of  intersection  nf  these  two  planar  curves  will  yield,  thrcigh  use  of 
(53),  all  possible  solutions  of  the  origmal  system  (52).  As  in  the  simpler  cases 
treated  previously,  the  nature  of  these  solutions  vill  depend  on  the  form  of  the 
functions  f2  as  well  as  on  the  side  conditions.  For  instance,  if  both  fj^  and 
f2  are  quadratic  forms,  then  (55)  reduces  to  a  single  equation  (in  either 
or  Xq)  of  the  fourth  order.  Hence,  if  no  pathological  cases  occur,  there  are, 
at  most,  four  distinct  solutions  to  the  original  problem. 

4.  The  preceding  discussion  points  the  way  to  a  more  general  extension 
of  the  same  type.  The  system  to  be  considered  now  is  of  the  form 

~  J  ^  -  j  '  .1  >  (56) 
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In  which  and  are  complex  numbers  and  the  vectors  are  still  assumed 
to  be  linearly  independent.  Use  of  the  linear  equp.cions  in  (56)  to  eli  .oinate,  say, 
the  first  n-r  Xj  values  reduces  system  (56)  to  a  system  of  r  equations 

)r,  (57) 

The  points  of  Intersection  of  these  r  hypersurfaces  will  then  yield  the  possible 
solutions  of  (56).  We  remark  that  ccnticuous  ma.tifolds  of  solutions  will  be  ob¬ 
tained  if  the  functions  are  not  linearly  independent. 

For  the  above  extensions,  one  can  formulate  corresponding  problems 

where  extrema  of  functions,  subject  to  linear  and  nonlinear  constxaints,  are 

required.  Suppose,  fur  example,  that  instead  of  system  (55),  we  must  determine 

the  extrema  of  FWx  ,,  x  )  subject  to  the  single  constraint 
n-i  n 


(58) 


This  can  be  done  either  by  straightforward  elimination  of  x^  ^  from  (57)  or  by 
the  use  of  Lagrange  multipliers. 


Linear  Dependence  of  Vectors  am 


We  shall  deal  now  with  the  extension  of  the  previous  results  to  the  case  of 
linear  dependence  of  the  vectors  a^  appearing  in  system  (1).  The  linear  de¬ 
pendence  of  these  vectors  implies  that  one  or  more  of  the  linear  equntioiis  In 
system  (1)  are  expressible  as  linear  combinations  of  the  rest.  If  this  is  the 
case,  the  redundant  equations  can  be  deleted  and  the  remaining,  say,  s  equa¬ 
tions  will  be  linearly  independent.  We  jhall  immediately  formulate  the  problem 
so  as  to  include  the  more  general  system 

X  —  m  ^ .  S  ^  S  -C  n-  r  .  (59) 
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As  compared  with  system  (56),  the  present  system  (59)  consists  of  fewer  equa¬ 
tions;  i.  e. ,  the  number  of  unknowns  exceeds  the  number  of  equations.  The 
vectors  ajQ  appearing  in  (5S)  arc  again  linearly  independent,  and  by  a  nimilar 
procedure  the  first  s  xj  values  can  be  eliminated.  The  resulting  system  will 
then  be 


The  original  assiunptioa  s  <  n-r  implies  also  that  r  <  n-s.  Therefore,  the 
possible  solutions  of  (60)  ’.vlli  depend  on  n-s-r  parameters,  and  the  solution 
of  the  original  system  (59)  may  include  continuous  manifolds  of  n-s-r  dimen¬ 
sions. 


The  simplest  system  that  can  be  solved  explicitly  is 

'px  -hX  Qx  — 

X  y  hn=  iy..  n-Z, 


(61) 


where  p  is  as  defined  in  system  (39). 

Assuming  that  we  already  have  the  correct  labeling,  we  can  write 

jr'f  =  ( i*  c*-x„  c*)/K\ 

where  all  the  entities  are  as  defined  in  (54). 

By  substitution  of  (62)  and  its  transpose  into  the  first  equation  of  (61), 
we  obtain 
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f */r’f c*  -K„Z’) ^ 

+ ( ^*-  Xr^  c/-x„  i3*A^'(  S'- 

-h  (L‘"- 

-t  =  h  . 


(63) 


All  the  starred  quantities  and  A  are  n-2-dimensional,  with  definitions  similar 
to  those  exhibited  in  (54),  while  q*,  q*  and  are  deClned  by 

X  ^ 


The  form  of  (63)  reflects  both  the  special  partitioning  of  the  variables  into  x* 
and  x(^)  and  the  assumed  symmetry  of  the  matrix  Q  of  the  quadratic  form  in 
system  (61). 

Since  Q*  is  also  symmetric,  we  can  write  (63)  as 
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where  the  coefficients  are  given  by 

in^fH  V  -1  '  '  i 

f>  =  C*^'  Q'^A-’^  +  yAA''^r  -r- l'^A~‘'c^ - 
t=  -  Zc^'^A^'c^-r-  c-A-'OTA-'c* 

^  -ZrA''^  -  zl^^'GTA-'^ 

T-  =  -  r'A'i^  -Z‘^A'^e -zi^^Q'^/r'cf 

^  =  ^/i-'^  -,  t*/r*QV|-'S  - 1. . 

(65) 

If  we  solve  for  from  (64)  we  will  get,  at  most,  two  distinct  solutions 

which  depend  parametrically  on  Xq.  The  original  problem  has  therefore,  at 
most,  two  distinct  one-parametar  manifolds  of  solutions  (if  no  pathological 
cases  arise). 

The  above  results  show  us  how  to  proceed  in  finding  the  extrema  of  the 
function 

subject  to  the  constraints 

^  ~  y>rr^  >  tn  =  '1^  ■  -  '  y  n  -2  . 

This  problem  is  equivalent  to  finding  the  extrema  of  the  function 

F (Xn^  jXn)—  -2P>Xn^Xr,  <68) 
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where  is  equal  to  the  (i  defined  in  (65)  with  b  =  0. 


If  we  assume  that  no  pathological  cases  arise,  there  is  a  single  extremum 
at  values  of  obtained  as  a  solution  to  the  system  of  the  following  linear 

equations : 

-r  C5-  —o 


namely. 


C 


r 


= 


s  -Zpx„^  : 

-h  ft  t: 


(69) 


^  ~  Z(  a  f) 


;  p>^-  o/  y  ^  o  . 


= 


(70) 


Substitution  of  these  values  into  (62)  gives  the  value  taken  by  the  vector  x*, 
and  *hus  the  solution  of  the  original  problem  has  been  found. 


V.  CONCLUSION 


In  the  precedin’  sections  we  have  presented  exact  algorithms  of  solution 
for  algebraic  systems  consisring  of  mixed  linear  and  nonlinear  equations.  These 
algorithms  involve  basically  the  evaluation  of  vaiious  determinants  as  well  as  a 
finite  number  of  products  and  summations.  This  in  itself  does  not  constitute  an 
obstacle  to  implementation  because  of  the  widespread  availability  of  computing 
facilities.  One  queutlon  that  can  be  raised  in  this  context  is  whether  these  algo¬ 
rithms  are  the  most  rapid  possible.  The  author  cannot  answer  this  equestlon 
since  alternative  exact  treatments  of  these  same  problems  do  not  seem  to  exist, 
and  therefore  there  is  no  basis  for  comparison. 

At  this  ooint  a  second  question  must  be  raised.  When  discussing  the  unique¬ 
ness  of  the  various  solutiv  ns,  we  have  assumed  that  a  nonzero  minor  ct^uld  be  found, 
f^ince  such  a  minor  is  not  given  from  the  start,  we  must  first  fint^  it  before  applying 
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the  algorithm.  For  example,  to  find  srch  a  minor  for  the  system  (59),  we  must 
examine  at  most  (”)  minors.  Obviously,  the  actual  number  exammed  may  be  con¬ 
siderably  smaller,  especially  if  several  such  aouzero  minors  exist.  This  brings 
into  prominence  the  second  question:  Does  the  application  of  the  algorithm  lead 
to  identical  solutions  when  various  nonzero  minors  are  used?  By  "identical”  we 
mean  that  a  change  in  labeling  will  bring  into  coincidence  solutions  based  on 
different  minors. 

The  answer  to  this  seems  to  be  in  the  affirmative;  i,  e.  ,  no  new  solutions 
are  Introduced  by  using  different  nonzero  minors  to  start  the  algorithm.  This 
can  be  verified  directly  for  systems  with  few  unknowns,  but  a  proof  for  a  general 
system  (in  particular  for  a  system  consisting  of  at  ’e;ist  one  highly  nonlinear 
equation)  is  not  known  to  the  author. 

Finally,  we  point  out  that  the  above  results  should  be  relevant  in  the  con¬ 
text  of  nonlinear  programming,'*  where  one  deals  with  systems  of  the  type  (59), 
ordinarily  with  r  =  1.  The  variables  x  =  (x^,  .  .  .,  Xq)  include  the  so-called 
slack  variables,  which  arise  from  the  transformation  of  a  system  of  inequalities 
to  a  system  of  equations.  The  given  single  function  f(X2,  ....  Xjj)  must  be  ex- 
tremized  subject  to  the  conditions  ajn'x^=  bj^^;  m  -  1,  .  .  . ,  s;  and  s  ;n-l.  Al¬ 
though  a  variety  of  methods  of  attacking  this  problem  have  been  in  use  for  years, 
it  is  important  to  realize  that  these  methods  yield  only  approximate  solutions. 
By  contrast,  the  direct  approach  described  here  will  always  yield  an  exact 
function  of  n-s-1  variables  to  be  extremized.  If  the  points  of  maxima  and  minima 
of  this  function  can  be  found  in  the  entire  n-s-i-dimensional  space,  ihen  the 
programming  problem  will  be  complecely  solved  if  use  is  made  of  the  residual 
coiLstralnts  of  non-negativity  of  the  x^^,  where  k  -  1,  .  .  .  ,  n. 

It  should  be  clear  from  this  discussion  that,  for  these  nonlinear  programming 
problems,  the  actual  choice  of  a  method  will  be  dictated  by  the  capacity  of  the  com¬ 
puting  facility  available,  the  computing  time  required,  and  the  penalties  involved 
in  the  use  of  approximate  extremal  points. 
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